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Abstract. In this short paper we present the main features of 
a new quantum programming language proposed recently by Peter 
Selinger which gives a good idea about the difficulties of constructing 
a scalable quantum computer. We show how some of these difficul- 
ties are related to the contextuality of quantum observables and to 
the abstract and statistical character of quantun theory (QT). We 
discuss also, in some detail, the statistical interpretation (SI) of QT 
and the contextuality of observables indicating the importance of 
these concepts for the whole domain of quantum information. 
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1.1 Introduction 

In classical information theory a fundamental unit is a bit taking the values 
or 1 which can be easily implemented by presence or absence of a particular 
physical signal. In quantum information theory a fundamental unit is a qubit 
whose values are vectors in a Hilbert space Hi=C 2 which are in one-to-one 
correspondence with the states of some two level physical system. Using two bits 
we may send 4 messages : 00, 01, 10, 11. With two qubits the messages are the 
vectors in H2=C 2 £g)C 2 with n qubits they are the vectors in H„. This richeness of 
the information together with long range correlations of the entangled quantum 
states stimilatcd a vigorous research on the quantum information theory and in 
particular on the quantum computer project. [10, 15, 30] 

The quantum computation has been usually studied in terms of the Turing 
machines , or in terms of circuits and gates. Several detailed quantum program- 
ming languages have been recently developped[7, 23, 24, 31, 38] . In particular 
an attractive functional programming language was proposed by Selinger [35] . 
With all these developments understand better what kind of the quantum hard- 
ware is needed . This quantum device (computer) has to manipulate the states 
of the qubits in a well defined rapid and precise way. These states have to 
to be implemented as the states of a particular quantum system. There are 
several proposed physical implementations of qubits and quantum gates using: 



trapped ions[12, 22] , cavity QED[40] , nuclear spins[21], Josephson junctions [32] 
and quantum dots[17, 29, 41] . The states of interest of the quantum register are 
all pure states of n qubits. To make a quantum computer, with error correcting 
ability, the quantum register should contain from 10 4 to 10 6 qubits . One should 
be able to to perform any unitary operation on the qubits using an appropri- 
ate sequences of some universal set of quantum gates [14,30,41] . One should 
be able to make a measurement in some chosen common measurement basis 
on each chosen qubit without disturbing the neighbouring qubits. The interac- 
tions of the system with the environment should be controlled in order to avoid 
a decoherence of the created pure states[4, 5, 20]. Finally the system should 
be scalable it means that the difficulty of performing gates , measurements , 
etc., should not grow too fast with the number of qubits in the registerin this 
short. In spite of the magnificent technological progress in nanophysics[17, 22] 
and optics [10, 28, 42]we are still far from the engineering stage of the quantum 
computer and the goal of 10~ 5 errors in the functioning of the gates is far from 
being achieved experimentally[39, 20] . It is relatively easy task, on a paper, to 
maniputate entangled states of the qubits, to switch on and off the interactions, 
to apply a sequence of the unitary operators or to perform a particular mea- 
surement leading to the instantcnuous collapse of the quantum state. To do it 
meaningfully and in a controllable way in the experiment is much more difficult 
because the quantum states are not the attributes of the individual quantum 
systems. 

The states of the classical objects are characterized by well defined attributes 
which may be changed in the deterministic way quasi instantaneously . The 
states of the quantum system are described by the vectors $ or by a density 
operator g which give the statistical predictions for the repeated measurements 
of quantum observables performed after the identical preparatory manipulations 
on the same or on the "identical" quantum systems. [8, 9, 16] . 

The problem is that our knowledge about the microworld is always indirect 
[11,25] and that we may only estimate the values of the physical observables 
from the empirical distributions of the experimental results. No experimen- 
tal result in the microworld may be predicted with the certainty [9, 27]and the 
quantum theory supplemented sometimes by some stochastic ad hoc assump- 
tions gives only the predictions on the statistical regularities observed in the 
experiments. This applies as well toall classic quantum experiments as to the 
experiments studying the interference effects using low intensity sources of en- 
tangled photons [33,34,42] or beams of heavy Ceo molecules[6]. This applies 
also to the experiments with the trapped ions and the quantum dots. 

It seems to us that any successful device in quantum information must take 
into consideration the fact that the values of all quantum observables are con- 
textual and the manipulation of the quantum state has nothing to do with the 
change of the attributes of some well defined and localised physical microsystem. 

The paper is organized as follows. In the section 1 we give a short discussion 
of the attributive and contextual observables in physics. In the section 2 we 
recall a formalism and the statistical interpretation of the quantum theory.In 
the section 3 we talk about qubits and quantum gates. In the section 4 we 
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present a principle features of the Selinger's quantum programming language. In 
section 5 we mention few implementation problems. 

1.2 Contextual observables 

Everyday observations taught us that the classical objects can be described and 
classified according to their characteristic properties, called attributes, such as 
their: form, size, dimension, colour etc. We found the effective classifications of 
plants, animals, chemical substances, stars and elementary particles. The static 
attributes are not sufficient to differentiate between different behaviour of the 
physical objects in various experiments this is why the contextual properties 
have been introduced. 

An attributive property as a constant property of an object which can be 
observed and measured at any time and which is not modified by the measure- 
ment. 

A contextual property as a property revealed only in a specific experiment 
and characterizing the interaction of the object with the measuring apparatus. 
Let us quote Accardi's example of a colour of a chameleon which is green on 
the leaf and brown on the bark of the tree what ressembles the behaviour of 
quantum systems [1,2]. 

Another important example of contextual properties , inspired by some 
random experiments we discussed in [26], is the following. Let us imagine that 
we have several double sided coins Ci,..Cfc having two faces : red ("R") and 
green ("G") . The coin my be flipped by different flipping devices Di...D„. 
Each time if we flip a coin Cj using the device D., we assume that we have an 
independent random Bernoulli trial with a probability of a success: p(R)=pij. 
If pij . T^Pfcz if i^k and j^l we clearly see that the probability py is neither an 
attribute of the coin Cmor the attribute of the device Dj . Each value pij may 
be called a contextual property of the coin Q. In general the values p,j are not 
known in advance and they can only be estimated within the error bars from 
the the long runs of the corresponding random experiments. 

A much more detailed discussion of the contextual observables may be found 
in [25] where we considered the sources of some hypothetical particle beams, 
detectors( counters), filters, transmitters and instruments and we obtained the 
following general results.: 

1) Properties of the beams depend on the properties of the devices and 
vice-versa and are defined only in terms of the observed interactions between 
them. For example a beam b is characterized by the statistical distribution of 
outcomes obtained by passing by all the devices d». A device d is defined by 
the statistical distributions of the results it produces for all available beams bj . 
All observables are contextual and physical phenomena observed depend on the 
richness of the beams and of the devices. 

2) In different runs of the experiments we observe the beams bk each char- 
acterized by its empirical probability distribution. Only if an ensemble fi of 
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all these beams is a pure ensemble of pure beams we can associate the esti- 
mated probability distributions of the results with the beams befi and with the 
individual particles members of these beams. 

3) A pure ensemble fi of pure beams b is characterized by such probability 
distribution s(r) which remains approximately unchanged: 

(i) for the new ensembles fi, obtained from the ensemble fi by the application 
of the i-th intensity reduction procedure on each beam befi 

(ii) for all rich sub-ensembles of fi chosen in a random way 

The quantum theory gives the probabilistic predictions thus it provides only 
a contextual description of the physical systems. 

There is a wholeness in any physical experiment [8, 9] .The initial states are 
prepared and calibrated, they interact with the experimental arrangement and 
the modified final states and/or the final numerical results the measurements are 
found. The quantum theory (QT) does not give any intuitive spatio-temporal 
picture of what is physically happening ,the QT gives only the predictions about 
the final states and about statistical distribution of the counts of the detectors. 

1.3 Quantum formalism. 

Let us recall now a standard description of the pure state of a quantum sys- 
tem with infinite number of degrees of freedom We consider a Gelfand triplet 
of spaces C H C Q! where is for example a Schwartz space of rapidly 
decreasing and infinitely differentiable functions on R n , H=L 2 (R n ) and Q! is a 
space of tempered distributions on ft. Let us consider two observables A and B 
measured in the mutually exclusive experiments I and II which are represented 
by self adjoint non commuting operators A and B. If {p^ } AeA and {<p^}^ gr 
are complete sets of the generalised cigenfunctions ( tempered distributions) of 
A and B then we get two different eigenvalue expansions of the unit state 
vector <I> e H of a studied physical system: 

*=J> A (A)^dA (1) 

A 

and 

*=J> B (7)<d 7 (2) 
r 

where ip A (^) — (<p£,ty) and V B (7) = (Pj^) are square integrable com- 
plex value functions. The probabilities that the measured value of A in the 
experiment I will fall into the interval [a,b] and that the measured value of B 
in the experiment II will fall in the interval [c,d] are given by the following well 
known formulae: 

P(a<A<b) = / A^jA) V, (A)dA , (3) 

a 
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P(c<B<d) = / 7 ^ B (7) ^ B (7)d7 • (4) 

c 

If a state of the system changes in time = U(t)^, where U(t) is a 

unitary time evolution operator. The formulae (1) and (2) are mathematically 
equivalent but they describe completely different experiments. The predictions 
of QT for any experiment measuring the values of the observable O are always 
obtained using a couple ( a state <J/ and an operator O ) and clearly all the 
information about the physical system is contextual and probabilistic. A value of 
a physical observable is not an attribute of the system but it is the characteristic 
of the pure ensemble of the experimental results created in the interaction of 
the system with the experimental device. 

The fact that using the same state vector ^ we may obtain the predictions 
on the results of all different measurements performed on the system led to 
the statements that the QT provides a complete description of each individual 
physical system. From this statement there is a one step to treat a state vector 
as an attribute of each single individual system what is incorrect and what 
leads to false paradoxes [9, 16, 27]. 

The most states in physics are not pure but mixed. In the case of a m- 
level system the most suitable is the formalism of the density operator /? acting 
in m dimensional hilbert space H m such that the expectation value of of any 
observable O is given by 

(0)=Tr(pO) . ^ (5) 
The time evolution of p(t) is given by: 

p(t) = u(t)pu(ty (6) 

where * denotes the hermitian conjugation. 

The density operator formalism is more general since by coupling a system 
to the environment one may describe in terms of the reduced dynamics [3, 19] 
the passage from the pure to the mixed states allowing a description of the 
decoherence or of the measurement process. One gets simpler formalism by 
representing the operator pbymxm positive hermitian matrix p. The density 
matrices and their transformations are building blocks of the quantum program- 
ming language [35]. In the follownig two sections we present the elements of this 
language. 

1.4 Qubits and quantum gates 

A pure state of a single qubit is a formal linear combination of two known states 
of some computational basis:q=aO + j31 which may be represented by a unit 
column vector u in C 2 and /or by 2x2 complex density matrix u u* . The 
n qubit states are spanned by the tensor products of n single qubit states and 
all their states pure and mixed may be represented by 2"x2" complex density 
matrices M.. 
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A unitary operations S on n qubit states are called n-ary quantum gates and 
their action on the state is defined by equation (6): SMS*. 

If S is n-ary quantum gate, then the corresponding controlled gate S c is the 
n+1- ary gate defined by the 2™ +1 x2" +1 matrix:: 



S c = 



I 

s 



(7) 



where I and are 2™ x 2™identity matrix and zero matrix respectively . In 
general a matrix with square brackets will be a matrix composed of 4 blocks 
each containing square matrices of the same dimensions. 

The following set of four unary and 5 binary gates can be chosen to be built 
into the hardware [35]: 



N= 



1 

1 



H =7i 



1 1 
1 -1 



V= 



1 
i 



1 

Vi 



and X= 



10 
10 
10 
1 



(8) 



A set is not uniquc[30, 41] but it is sufficient to approximate the action of 
any unitary gate according to the theorem [14]which says that for any unitary 
matrix 

Se C 2 " x2 " and e ) there exist a unitary matrix S' and a unit complex 
number A such that 

||S-AS'||<e (9) 



and such S'=I®A®J where I, J are identity matrices of the appropriate di- 
mensions, and A is one of the gates H, V c or X. 

In supplement to the unitary reversible operations, the operation " measure" 
is introduced: 
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(10) 



where the first matrix corresponding to a classical bit is obtained with the 
probability Tr A and the second corresponding to the classical bit 1 is obtained 
with the probability Tr B . If the classical bits of information are ignored the 
final result is a normal sum of the outcome matrices(lO): 



" A 


B 




' A 





C 


D 







D 



(11) 



G 



1.5 Quantum programming languages 



Several proposals for the quantum programming languages has been made. Let 
us mention here the pseudo-cod formulation by Knill [23,24], Omer's rich pro- 
cedural QCL[31], Bottom's quantum C++ [7]and Sanders and Zuliani qGCL 
[38] .All these quantum programming languages are so called imperative type 
languages because the quantum data are manipulated in terms of arrays of quan- 
tum bit references which requires the insertion of a number of run-time checks 
into the compiled code which must be executed at run- time not at the compile 
time. 

The quantum programming language proposed by Selinger[35] is a functional 
programming language with a static type system which guarantees the absence 
of any run-time errors. The syntax and semantics of the language allow high- 
level features such as loops, recursive procedures, and structured data types. 
Each statement operates by transforming a specific set of inputs to outputs, 
and the principle of non-cloning of quantum data is enforced by the syntax. Both 
data flow and control flow are described using the paradigm "quantum data, 
classical control" in consistence with Shor's factoring algorithm, Grover's search 
algorithm and the Quantum Fourier Transform [13, 14, 30, 36, 37]. 

The models contains the classical and quantum flow charts . Quantum flow 
charts are similar to classical flow charts, except that a new type qbit of quan- 
tum bits, and two new operations:unitary transformations and measurements 
are added.. The following notation is used for these new operations: q*=S for 
application of a unary unitary quantum gate S to a quantum bit q, p,q*=S for 
a binary quantum gate S applied to a pair p ,q of quantum bits, < measure p>: 
for the branching statement giving a couple of outcomes after the measurement 
of p, "o" for the merge operation yielding the final mixed state etc. 

In an an example taken from [35] a flow chart below corresponds to the 
fragment of the program which: inputs two quantum bits p and q, measures 
p, and then performs one of two possible unitary transformations depending on 
the outcome of the measurement. The output is the modified pair p,q. 

input p,q: qbit 
I 

measure p 
0/ \1 

q*=N p*=N (11) 

\ / 
merge (o) 

I 

output p,q: qbit 

If the input is a general mixed state described by a 4x4 density matrix M 
written in the block notation the flow chart describes the following sequence of 
operations performed on the density matrix M. 
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where A' = NAN*+D . Let us note that by the usual mormalisation con- 
vention the traces of all the density matrices are equal to the probabilities that 
the corresponding edge is reached. If the input is a pure state, then the states 
along each of the two branches of the measurement continue to be pure and a 
final mixed state is obtained by the merge operation.. 

1.6 Hardware requirements 

The quantum computing language could be implemented on various devices but 
a real computation time gain could be only achieved using a quantum computer 
par example a QRAM type machine [23] which consists of a classical computer 
which controls a special quantum hardware device such that [35]: 

a) The quantum device provides a potentially large number of individually 
addressable quantum bits. 

b) Quantum bits can be manipulated via two fundamental operations: uni- 
tary transformations and measurements. 

c) The quantum device will implement a fixed, finite set of unitary transfor- 
mations which operate on one or two quantum bits at a time. 

d) The classical controller communicates with the quantum device by sending 
a sequence of instructions, specifying which fundamental operations are to be 
performed. 

e) The only output from the quantum device consists of the results of mea- 
surements, which are sent back to the classical controller. 

The operating system has to keep a list of quantum bits that are currently 
in use by each process. When a process requests a new qubit, the operating 
system finds a qubit that is not currently in use, marks that qubit as being in 
use by the process, initializes the contents to 0, and returns the addresse of the 
newly allocated qubit. The process can then manipulate the qubit, for instance 
via operating system calls which take the qubit's addressas as a parameter. The 
operating system ensures that processes cannot access qubits that are not cur- 
rently allocated to them - this is very similar to classical memory management. 
Finally, when a process is finished using a certain qubit, it may deallocate it via 
another operating system call; the operating system will then reset the qubit 
to and mark it as unused. In practice, there are many ways of making this 
scheme more efficient, for instance by dividing the available qubits into regions, 
and allocating and deallocating them in blocks, rather than individually. 

Reseting or initializing a qubit to is implemented by first measuring the 
qubit, and then by performing a conditional "not" operation N c dependent on 
the outcome of the measurement. 
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In the introduction we reviewed various difficulties in constructing a scalable 
quantum computing device. Let us add here some additional problems related 
to the contextuality and SI of QT. 

1) Concerning the measurements (1 0).Let us consider an entangled multiqubit 
state* . The QT gives us statistical predictions on the results of various long 
runs of coincidence type measurements which may be performed by distant 
detectors on various qubits of a system.. For any particular qubit we may 
obtain or 1 with a given probability. Imagine now that in a given moment 
of time a single measurement performed on the qubit q gives the result 1 then 
according to SI we do not have any deterministic information about the results 
of other coincidence measurements performed in the same or future time on 
other qubits[9, 27]. The collapsed state vector gives us only the conditional 
probabilities for the the results of the coincidence measurements performed on 
the remaining qubits corresponding to subcnscmble of all n qubit events in which 
a measurement performed on the qubit q gave 1. 

2) Concerning the unitary gates (8,9). These transformations have to be 
implemented by sequences of: qubit -external field and/or qubit-qubit physical 
interactions. The ideal model requires an instataneous switching on and off the 
various physical interactions without disturbing the other qubits In practice one 
wants to use so called adiabatic switching which seems to be very difficult task. 

3) The produced pure states of quantum register have to be protected against 
the decoherence due to the incontrollable influence of the environment. 

1.7 Conclusions. 

In our opinion any successful design of a quantum hardware for the purpose 
of quantum information must be consistent with the contextual and statistical 
character of quantum theory. 
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